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When a single magnetic adatom is deposited on a surface of a metal, it affects the charge and spin 
texture of the electron gas surrounding it. The screening of the local moment by conduction electrons 
gives rise to the Kondo effect. Here we investigate the effect of the Rashba spin orbit coupling on 
the local magnetization density of states (LMDOS) around a Cobalt impurity on a Au(lll) surface 
in a magnetic field. We show that the in-plane component of the LMDOS is exclusively associated 
with the Rashba spin orbit interaction. This observation can be experimentally exploited to confirm 
the presence of the Rashba effect on surfaces, such as Au(lll), by performing spin-polarized STM 
measurements around the Kondo impurity. 

PACS numbers: 72.15.Qm, 72.25.-b, 75.70.Tj 
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I. INTRODUCTION 

In two-dimensional structures, the coupling between 
the spin and angular momentum can lead to a vari- 
ety of interesting phenomena.^ Interestingly enough, the 
Rashba spin-orbit (SO) interaction^ has become one of 
the intriguing and desired ingredients of modern na- 
noelectronics and spintronics.^ Studying the Rashba- 
induced effects for atoms placed on surfaces is especially 
interesting, because on one hand, it opens up a way 
to probe the strength of the SO coupling, and on the 
other hand, it may allow one to explore the interplay 
between the Kondo screening and the SO interaction. 
In this regard, the Au(lll) surface is very promising, 
since it exhibits a measurable energy splitting of sur- 
face band states, which was first experimentally observed 
by LaShell et al. [4]. However, because the Friedel os- 
cillations induced by an adatom remain practically un- 
affected by the Rashba SO interaction^, the presence 
of the Rashba effect could not be verified by typical 
scanning tunneling microscopy (STM) measurements, fo- 
cused mainly on the charge transport spectroscopy. Con- 
sequently, not much attention has been paid to it so far. 
It has been noticed only very recently that placing a mag- 
netic adatom on the surface can help probe the Rashba 
spin-orbit effect. By treating the impurity at the classi- 
cal level and in the absence of the Kondo effect, Lounis 
et al. [6] have shown that the induced spin polarization 
of the electron gas surrounding the magnetic adatom ex- 
hibits a spin texture, which is a superposition of two 
skyrmionic waves with opposite chirality. This has been 
attributed to the presence of the Rashba SO interaction. 

In the present work we pursue this problem further 
and study the local texture of the spin-resolved density 
of states around the magnetic adatom in the case of non- 
trivial many-body interactions, such as the ones leading 
to the Kondo effect,^ and in the presence of an exter- 
nal magnetic field applied along the zth direction. The 
Kondo effect, which occurs when a local impurity spin in 
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FIG. 1. (Color online) Spatial distribution of the radial 
component Mr (left panel) and of the total LMDOS M (right 
panel) at energy E — —2 meV below the Fermi surface, and 
for a magnetic field 5 = 3 T applied along the zlh direction. 
The smaller (orange) spheres indicate the position of the Au 
atoms on the Au(lll) surface, and the large (black) sphere 
represents the magnetic adatom. The arrows indicate the 
amplitude and direction of the LMDOS. 



a metallic host is screened by the conduction electrons, is 
undoubtedly one of the fundamental effects in condensed 
matter physics.^ Here we investigate how the induced lo- 
cal magnetization density of states (LMDOS) around a 
magnetic adatom in the Kondo regime, in the presence 
of an external magnetic field, is affected by the Rashba 
SO interaction. 

In the many body formalism, the LMDOS, M(r, £^), 
can be expressed in terms of the retarded, single particle 
Green's function, ^(r, r',£^), as: 



M(r,£;) 



-^mTr,{aa(r,r,^)}, (1) 

TT 



where cr are the Pauli matrices, E is the energy mea- 
sured with respect to the Fermi energy £^i?, and r is 
the in-plane distance from the impurity. The real-space 
Green's function ^(r, r,£^) shall be computed in terms 
of the many body T-matrix for the conduction electrons, 
which describes the scattering of the surface electrons off 
the impurity. 



To make quantitative estimates, we focus on a Co atom 
on a Au(lll) surface,^' ^^ for which a Kondo temperature 
Tk of about 75 K was extracted from STM spectroscopy 
measurements.^^ The Au(lll) surface is modeled within 
the tight binding approximation (TBA), which, in spite 
of its simphcity, is able to properly describe the disper- 
sion of the Au surface states. ^^ To capture the Kondo 
physics correctly, the Co impurity is described in terms 
of the Anderson model. ^^ The many body T-matrix is 
related to the Green's function describing the local or- 
bitals of the Co ion, (see Sec. II B 2) which can be then 
computed with the aid of the numerical renormalization 
group (NRG) method, ^^ known as the most versatile and 
accurate in treating quantum impurity problems. More- 
over, to make realistic predictions, in NRG calculations 
we take into account the full energy dependence of the 
density of states (DOS) of the Au(lll) surface. While in 
general the magnetic impurity itself can have a compli- 
cated orbital structure, and channels with different sym- 
metries may couple to the surface, within the NRG ap- 
proach the coupling is assumed to have an 5-wave sym- 
metry. 



II. THEORETICAL FRAMEWORK 

A. Modeling the Au(lll) surface 

Let us introduce the details of the lattice under in- 
vestigation, the tight-binding Hamiltonian describing it, 
and the corresponding band structure. The Au(lll) 
surface presents a hexagonal structure, with one atom 
per unit cell. The basis vectors of the direct lattice are 
ti = (\/3/2, 1/2) a and t2 = (-\/3/2, 1/2) a, with a the 
lattice constant (a = 5.75 A). Here we are particularly 
interested in the changes induced locally by a magnetic 
impurity in the LMDOS. We shall not address the so- 
called herringbone reconstruction^^, which may be rel- 
evant when analyzing photoemission spectra. Also, the 
external magnetic field is assumed to produce no kinetic 
effects on the surface states, as its effect is marginal. In 
spite of its simplicity, this tight binding description is 
rather robust, and can be checked against more sophis- 
ticated ab-initio band structure calculations,^^ or com- 
pared to experimentally measured binding energies.^ 



One of the main results of this paper-a non- vanishing 
in-plane magnetization Mp- is sketched in Fig. 1 together 
with the total magnetization M in a magnetic field, B = 
3 T. Notice that at 5 = 0, the local polarization vanishes, 
M = 0. Although the total LMDOS M = (Ma,,M^,M^) 
is an important quantity, we have found that the in-plane 
component Mr is much more interesting, as it is strongly 
affected by the presence of the SO interaction. On the 
other hand, the out of plane part, M^, depends weakly 
on the Rashba interaction and displays a spatial behavior 
that is somewhat similar to the one observed in the local 
density of states (LDOS),^'^"^'^^ More than that, Mr is a 
pure Rashba effect, as it vanishes if the SO interaction is 
turned off. 



Experimentally, it is possible to measure the radial 
component of the energy-dependent LMDOS with the 
state-of-the-art spin-polarized STM techniques, ^^ which 
can thus provide an important information on the pres- 
ence and strength of the Rashba SO interaction. In this 
paper, using realistic parameters, we study the behav- 
ior of the energy and position-dependent LMDOS in the 
region around the magnetic impurity. 



The paper is organized as follows: In section II we 
introduce our model Hamiltonian that describes the Au 
surface. The description of the magnetic impurity prob- 
lem is also presented in the same section. In section III 
we present and discuss our numerical results. We close 
with conclusions in IV, where our main findings are reit- 
erated. 



1. Hamiltonian for the Au(lll) surface 

We model the Au(lll) surface in terms of a tight 
binding Hamiltonian,^^ taking into account the hopping 
between the nearest-neighbor p^-orbitals subject to the 
Rashba SO interaction 



Ho= Yl I^ (^r,r' + ^p (^r,rO <a^r 



(2) 



<r,r'> (7 

*" Y. Y. 

<r,r'> aa^ 
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^L^r',a' . 



The first term describes the hopping and the on-site en- 
ergies, while the second one is due to the Rashba spin- 
orbit coupling. Here, ^J ^ creates an electron in the Au 
p^-orbital at position r with spin a, try are the nearest- 
neighbor hoppings between these orbit als, and Ep de- 
notes their on-site energies. In the second term, a is the 
strength of Rashba interaction. 

By fitting the tight binding dispersion, £r(k), to the 
experimentally measured binding energy of Ref. [4] along 
the TM direction, one can extract the band parame- 
ters^^i t = -0.30 eV, Ep = 1.41 eV and a = 0.02 eV. 
We note that the effect of the external magnetic field on 
the surface electrons is rather minimal. A simple analysis 
of the energies involved, shows that the Zeeman splitting 
for a magnetic field of about 3 T (the g-factor was taken 
to be ^ = 2.5) is AEz = 0.43 meV, i.e. five times smaller 
than the Rashba energy: Er = -?>a^/{2t) = 2.08 meV, 
and tiny as compared to the band parameters t and Ep. 
Consequently, the effect of the magnetic field on the con- 
duction electrons is neglected. 

The Hamiltonian (2) can be diagonalized in Fourier 



space by expanding the field operators ^r,a as 



*r 
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Vn 



6 (^Tja v^j Ck,r 1 



(3) 



k,r 



where Q is the number of unit cehs, and Ck,r annihilates 
an electron with momentum k in the chiral band r = ±1. 
Then, the dispersion is: 6:^(k) = Ep ^ G{k) + r |F(k)|, 
with 



G(k) = 2t 2 cos -—kxa cos ( -/c^^a j + cos{kya) 
F(k) = -a (1 + V3i) sin I — A;^a + -kya j + 
+ (1 — v3z) sin 



2 2 ^ . 



2 sin kya 



(4) 

Here, the form factor G(k) is coming from the nearest- 
neighbor hopping and F(k) is due to the Rashba SO 
interaction. The chiral band energies £r(k), and the wave 
function amplitudes Cr^aO^) are determined by solving 
the eigenvalue equation 



^Haa' (k) er,a'(k) = SrO^) er,a(k) 



(5) 



H{k) 



(6) 



with iJo-(T'(k) the components of the Hamiltonian (2) in 
Fourier space and spin basis. Its matrix form is 

Ep + G(k) F(k) 
F(k)* Ep + G(k) 

The corresponding eigenvectors can be evaluated analyt- 
ically. 



(7) 
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e (k) - ^ (l -:^» 



Then, in terms of the creation and annihilation operators 
for electrons in the chiral basis, Hq becomes 



Ho = ^ ^r(k) cl^ci^,T, 



(8) 



with the operators Ck,r satisfying the canonical anti- 
commutation relations, {c\^^rjcl.f r'} ~ ^(^ ~ k')^r,r'- In 
this way, the surface can be described in terms of free 
states, but with some chiral band structure. 
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FIG. 2. The normalized DOS for the Au(lll) surface ob- 
tained within the TBA. The two insets show the singularities 
at the band tails due to the Rashba SO interaction. 



where 6{t) is the Heaviside function. Because of the spin- 
orbit interaction, Q^J,{ri^t;rj^t') is diagonal in the chi- 
ral, but not in the spin space. In this context, the Fourier 
transformed, non-interacting retarded Green's function 
for the chiral band r is: Qr (k, cj) = {uj — er{k) -\-iO~^)~^ . 
Then, transforming back to the spin space one gets 



GaM^i^^J^^)=Yl 



k,r 



e...(k)e;,,(k)e-'-(---- 
oj - SrCk) + iO+ 



(10) 



This expression allows us to evaluate the density of states 
(DOS) for the conduction electrons as: 



JO) 



1 



H = — Tr, ^ma^^j(r„r„a;). (11) 



In Fig. 2 we show the DOS computed by using Eq. (11) 
for the Au(lll) surface. It displays a large van Hove 
feature at cj ?^ 2.03 eV and two sharp singularities at the 
band tails, see the insets in Fig. 2. The latter features are 
induced by the presence of the Rashba SO interaction. 
This DOS will be the input for the NRG calculations 
when solving the quantum impurity problem. 



2. Non-interacting Green's function 

In the non- interacting limit, the retarded Green's func- 
tion in the real space is defined in terms of the field op- 
erators for the conduction electrons as, 

C;fj,(r„i;r„i') = -i0(i - i'){{*r„.(t), <,.,(t')}>o, 

(9) 



B. Modeling the quantum impurity 

1. Impurity Hamiltonian 

To carry out the quantitative analysis of our magnetic 
impurity problem, we first need to establish how the mag- 
netic ion couples to the chiral bands. We consider here 
the top configuration^ in which one Co atom is located 



on top of an Au atom, to which it hybridizes. The hy- 
bridization with ah the other neighboring Au atoms is 
neglected. We have also considered other geometrical 
configurations (results are not presented here), where for 
example the Co atom is placed in plane, in the middle of 
a hexagon, or substitutes an Au atom, hybridizing with 
the nearest neighbors. Despite its simplicity, the con- 
sidered top configuration captures entirely the essential 
physics. Correspondingly, the hybridization Hamiltonian 
is written as 



Hv = vY,(^^l^^^,d, + h.c.). 



(12) 



Here rimp labels the Au site below the Co impurity, do- 
annihilates an electron with spin a at the Co orbital, and 
V denotes the hopping between the two orbit als. Within 
the NRG approach, it is convenient to model the mag- 
netic impurity by a single local orbital which carries only 
a spin label. Transformed to the chiral basis, Eq. (12) 
becomes 



HV=J2Y1 {VrA^) Ck.r*^- + h.C. 



=± k,, 



(13) 



This expression is quite general, and the particular loca- 
tion of the impurity atom is reflected in the k-dependence 
of the hybridization factor, Vr^aO^)- For the top config- 
uration, one has K,cr(k) = Fe*^^(k). Using Eq. (7) for 
the eigenvectors, one then finds: ^/^-^^(k) = V/V2 and 
Vr,i{k) = r VF(k)/(v^|F(k)|). With the hybridization 
Hamiltonian (13) at hand, the total Hamiltonian that 
describes the Co ion itself and the hybridization to the 
surface is: 



H = Hin 



Hv 



with 



H^ 



imp 






£d( 



ydld^ ^ Un^ni 



(14) 



(15) 



This Hamiltonian is similar to the single-impurity An- 
derson Hamiltonian,^^ but with a somewhat modified hy- 
bridization. The first term in Eq. (15) describes the on- 
site energy Sda = £d ^ \gi~^BB of the localized orbital, 
where we included a Zeeman splitting term due to the 
external magnetic field B applied along the zih direc- 
tion. We assume that the g-factor for the Co atom on 
the Au(lll) surface is around ^ = 2.5. In the second 
term, U represents the Coulomb repulsion felt when two 
electrons with opposite spins occupy the orbital, with 
Ua- = dl^da- denoting the occupation number. We take U 
and Sd from ab-initio calculations^^: Sd = —0.84 eV and 
U = 2.85 eV. The hybridization amplitude V is fixed 
by the Kondo temperature itself. Here, we define the 
Kondo temperature Tx as the half width at half maxi- 
mum (HWHM) of the spectral function for the local or- 
bital operator d^ in the absence of an external magnetic 
field. Then, to get Tk = 75 K, we take V = 0.25 eV. 



In the presence of a 5'[/(2) spin symmetry, the electrons 
in the spin={t, i} channels are scattered in the same way 
by the magnetic impurity. In order to observe any spin- 
resolved signal, it is necessary to break this symmetry by 
applying an external magnetic field along the zth direc- 
tion. Then, the scattering becomes spin dependent, as 
the Kondo resonance is spin-split. ^^'^^ One drawback of 
such a set-up is due to the large Kondo temperature^^: 
a relatively large magnetic field is necessary in order to 
produce a detectable splitting of the Kondo resonance. 
Here we have considered B = 3T. 



2. Calculation of the T -matrix 

To solve the quantum impurity problem, we employ 
Wilson's NRG method. ^^ NRG is a powerful tool for 
accurate calculations of equilibrium properties of arbi- 
trarily complex quantum impurities coupled to electron 
reservoirs. ^^ The method consists in the logarithmic dis- 
cretization of the continuum of conduction states, fol- 
lowed by a mapping to a one dimensional chain Hamil- 
tonian (Wilson chain) with exponentially decaying hop- 
pings. The mapping starts with expanding the operators 
Ck,r in terms of the eigenfunctions of the angular momen- 
tum^^, 



1 



Ck,r 



E 



1 



27r 



imcf) m 



(16) 



and then constructing an effective impurity model by in- 
tegrating out the electronic angular momentum modes. 
The broadening felt by the impurity is given by the imag- 
inary part of the hybridization function 



A.H= EE 



IK,.(k)|^ 



£^(k) + iO+' 



(17) 



To a first approximation, the Rashba spin orbit coupling 
is weak and does not affect the impurity spectral func- 
tion: A^(co') = A|(c(j). 

In order to determine the LMDOS, Eq. (1), one needs 
to calculate the full Green's function ^crcr/(r; r;a;), which 
can be expressed in terms of the T-matrix by using the 
Dyson equation. Let us now focus on the calculation 
of the T-matrix itself. For quantum impurity models, 
one of the most elegant ways to perform this task is to 
relate it to some local correlation function that can be 
computed numerically with the NRG. For the Anderson 
model, the T-matrix is related to the Green's function of 
the da- operators. ^^'^^ For a constant and real coupling 
V, the imaginary and real parts of the spin-resolved T- 
matrix are then given by 



^m T(^) 
3?e %{ou) 



-nV^Aaiuj), 



V^V 



I 



dw 



,Aa{u}') 



(18) 



with A(j{oj) the spectral function of the d^ operators and 
V denoting a principal value integral. In order to obtain 
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FIG. 3. (Color online) The energy dependence of the spin- 
dependent spectral function Aa{(^) of the local orbital around 
the Fermi level, calculated within the NRG approach. The 
splitting of the Kondo resonance is due to the external mag- 
netic field, 5 = 3 T. The inset shows the full energy depen- 
dence of Aa{(jo). The parameters for the Anderson model are: 
V = 0.25 eV, U = 2.85 eV and Sd = -0.84 eV. 



reliable results for the spin-dependent spectral functions, 
we have employed the density-matrix NRG.^^ In addi- 
tion, we have included in our calculations the full energy 
dependence of g'^^^E) for the Au(lll) surface. 

In Fig. 3 we show the energy dependence of the spin- 
resolved spectral function Aa{oj) in the vicinity of the 
Fermi level, with the inset presenting its full energy de- 
pendence. Although g'^^^E) shows a gap below E < 
0.42 meV (the bandwidth in the NRG calculations was 
fixed to D = 2.5 eV), the two Hubbard satellites and the 
Kondo peak at the Fermi level are clearly visible. The 
splitting of the Kondo resonance for 5 = 3 T is visible in 
Aa{oj). The applied magnetic field is not strong enough 
to suppress the Kondo resonance, however it is sufficient 
to produce a spin- resolved response detectable in the LM- 
DOS. The Aa{^) in Fig. 3 was computed at T = 0, but 
it can be argued that our findings remain valid as long 
as we are in the Kondo regime: T < min{Tx,5}. If 
the temperature increases, the Kondo peak becomes sup- 
pressed and eventually, at high temperatures T :^ Tk, it 
is completely smeared out by thermal fiuct nations. 



III. LMDOS: ANALYSIS OF THE NUMERICAL 
RESULTS 

In this section we shall describe how we compute the 
LMDOS. As discussed in Sec. I, the LMDOS can be 
related to the single particle Green's function, see Eq. 
(1), and satisfies the Dyson equation, when expressed in 
terms of the T— matrix. While this expression is some- 
what cumbersome in the chiral basis due to the presence 







12 3 4 

X [lattice constants] 





-0.2 


-0.1 0.0 0.1 0.2 

Energy [eV] 




d 


impurity 


6 


L^ 




, 


1 




nearest neighbor 


S' 


A 


03 






2".1 


r 






next nearest neighbor 


-1 


'^f 



Energy [eV] 



FIG. 4. (Color online) (a) The spatial distribution of the 2;th 
component of the LMDOS, M^, around the Kondo impurity 
at energy E = —2 meV, below Ef and (c) the corresponding 
cuts at different energies. The impurity is located at the cen- 
ter. The pattern formed around the magnetic impurity has 
a hexagonal symmetry and exhibits oscillations with distance 
from the impurity, with a period depending on the energy. 
The bright (dark) color corresponds to maximum (minimum) 
value of M^. (b) M^ as a function of energy and distance 
from the impurity along the x direction, and (d) the energy 
dependence of M^ (E) calculated at the impurity site and the 
two nearest-neighbor sites. 



of different form factors, in the spin space it simplifies 
considerably: 

Gaa' (r, r,uj) = 5^^l (r, r', cj) + ^ ^^^l (r, rimp, cj) x 

Ta- (rimp, cj) 0^?^/ (rimp, r', cj). (19) 

The LMDOS can be calculated from (5 ^crcr' (r, r, cj) = 
^^^/(r,r,cj) - ^^^l(r,r,u;), by replacing g(r,r,cj) -^ 
(5^(r, r,cj) in Eq. (1). Notice that the spin impurity acts 
as a simple point scatterer, and that in the magnetic re- 
sponse, only the second term in Eq. (19) gives a finite 
contribution. 

Within the present formalism, we are able to analyze 
both the spatial and the energy dependence of the lo- 
cal magnetization components. Here we will focus on 
Mr and the out-of-plane, M^ component of the LMDOS. 
The spatial distribution of the M^, calculated at energy 
E = —2 meV, is displayed in Fig. 4(a). The bright (dark) 
areas correspond to the maxima (minima) of M;^. One 
can see that its spatial dependence displays a hexagonal 
symmetry with respect to the position of the magnetic 
impurity. Moreover, M^ exhibits oscillations with an en- 
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FIG. 5. (Color online) (a) The spatial distribution of the in- 
plane X component of the LMDOS, M^, around the Kondo 
impurity at energy E = —2 meV, and (c) the corresponding 
cuts (in the x direction) at different energies, (b) Mr as a 
function of energy and distance from the impurity along the x 
direction, and (d) the energy dependence of Mr(c(;) calculated 
at the nearest and next-nearest neighbor sites. 



ergy dependent period. This behavior is displayed in 
Fig. 4(c). When moving away from the Fermi energy, 
the magnitude of M^ decreases and the period of the 
oscillations becomes shorter as the energy is increased 
from negative to positive values. This is presented in 
Fig. 4(b), which explicitly shows the energy and spatial 
dependence of M^ when E is swept across the Fermi sur- 
face. Close to Ep^ the asymmetry induced by B in the 
spin sector, together with the presence of the split Kondo 
resonance, maximizes the amplitude of the local mag- 
netization. However, when the energy is detuned from 
Ef, the LMDOS becomes suddenly suppressed, leading 
to a Fano-like resonance, similar to the LDOS.^^ More- 
over, the shape and magnitude of such a Fano resonance 
changes as one moves away from the impurity site, see 
Fig. 4(d). One should note that within the TBA, al- 
though we have represented the spatial distributions as 
continuous, the calculations are only valid at the atomic 
sites. 

In a finite 5, the features observed in M^ are present 
irrespective of the presence of the Rashba SO interaction 
at the surface. On the other hand, M^ is very sensitive 
to the Rashba effect. Finite a implies M^ ^ 0, otherwise 
M^ vanishes. The spatial distribution of M^ around the 
magnetic impurity together with its energy dependence 
are presented in Fig. 5. First of all, one can note that 
in the vicinity of magnetic impurity, the amplitude of 



M^ is smaller by approximately one order of magnitude 
than the amplitude of M^. In fact, M^ vanishes exactly 
at the impurity site, while M^ has a maximum there, 
see for example Fig. 4(a,c) and Fig. 5(a,c). However, 
with increasing distance from the impurity, x > 3 ^ 4a, 
both M^ and M^ become comparable and, in fact, for 
larger distances, the radial component can overtake M^, 
as its decay is much slower than that of the zth com- 
ponent. The basic properties of the spatial dependence 
of Mr can be deduced from Fig. 5(a). One can see that 
Mr is an odd function with respect to the radial distance, 
and oscillates with approximately the same period as M^. 
Again, the highest amplitude of M^ occurs for energies 
around the Fermi energy due to the Kondo effect, and 
as the energy increases, the period of the oscillations de- 
creases, see Fig. 5(b) and Fig. 5(c). The spatial and 
energy dependence of M^ is shown in Fig. 5(b), while M^ 
at some particular positions is displayed in panel (d). It 
can be seen that the shape of the resonance at the nearest 
neighbor and the second-nearest neighbor sites is similar 
to that of M^ [cf. Fig. 4(d)]. 

Finally, we would like to emphasize that M^ is on one 
hand proportional to the strength of the SO interaction, 
and on the other hand to the asymmetry of the T-matrix 
components for the spin-^ and spin-| channels. While a 
is an intrinsic feature of the surface, the asymmetry be- 
tween the spin-t and spin-| channels can be changed by 
simply applying an external magnetic field. This guar- 
antees that the topographic map of the surface devel- 
ops interference patterns in M^ if the Rashba interaction 
is present. Therefore, the measurement of the in-plane 
component of the LMDOS offers an alternative way to 
angle resolved photoemission spectroscopy (ARPES)^ to 
identify surfaces with spin orbit interaction. 



IV. CONCLUDING REMARKS 

In the present work we have investigated the behav- 
ior of the local magnetization density of states around 
a magnetic impurity in the Kondo regime, coupled to a 
metallic surface with Rashba spin orbit interaction. In 
order to make realistic estimates, we have considered a 
Co impurity on top of a Au(lll) surface. This problem 
has been addressed using band structure calculations and 
the NRG method, which allowed us to obtain reliable pre- 
dictions for the LMDOS. 

In particular, we have studied the spatial and energy 
dependence of the LMDOS for both the radial (in-plane) 
and zth (out-of-plane) components. We have found that 
the in-plane component of the LMDOS is a pure Rashba 
effect. Furthermore, it turned out that the radial com- 
ponent displays oscillations with the distance from the 
impurity, and decays much slower than the zth compo- 
nent, so that at larger distances, the in-plane compo- 
nent may become dominant. Since M^ vanishes in the 
absence of spin orbit interaction, measuring the radial 
component using spin-polarized STM, provides a way to 



confirm/ infirm the presence of the Rashba effect on sur- 
faces. Our observations provide thus an alternative route 
to investigate the spin splitting of any surface states, first 
observed by angle resolved photoemission spectroscopy^ 
in Au(lll). 
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